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Abstract. We propose a method by which to examine all possible partial 
difference Lax pairs that consist of two 2x2 discrete linear problems, where 
the matrices contain one separable term in each entry. We thereby derive 
new, higher-order versions of the lattice sine-Gordon and lattice modified KdV 
equations, while showing that there can be no other partial difference equations 
associated with this type of Lax pair. 



1. Introduction 

Nonlinear integrable lattice equations provide a natural discrete extension of clas- 
sically integrable systems. An extensively studied example is the lattice modified 
Korteweg-de Vries equation, referred to as LMKdV: 

{pxi+i, m -rxi, m+1 ) 
LMKdV : x l+x , m+1 = xr m -, r- 

[PXLm+l - rXl + i tm ) 

where xi tTn is the dependent variable, I and to are the discrete independent variables, 
and p and r are parameters. This equation provides an integrable discrete version 
of the well-known modified Korteweg-de Vries equation. Note that, in the past, 
the parameters p and r have been taken to be constant pQ, although singularity 
confinement has provided a way to de-autonomize LMKdV while maintaining its 
integrability, by allowing p and r to depend on I and to in a specific way [Sj . The 
integrability of such equations lies in the fact that they can be solved through an 
associated linear problem called a Lax pair. 

Lax pairs can appear in many guises, the type that we are exclusively concerned 
with in this paper consist of a pair of linear problems written: 

d{l + l,m) = L(l,m)0(l,m) 

e(l,m+l) = M(l,m)9{l,m). [ ' 

where 9(1, to) is a two-component vector and L(l, to) and M(l, to) are 2x2 matrices. 
These linear problems are described by L and M, which are referred to as the Lax 
matrices. If we were examining an explicitly determined Lax pair for LMKdV for 
example, the variable xi >m , along with the parameters p and r, would be found to 
reside within L and M. The easily derived compatibility condition on this Lax pair 
is 

L(l,m + l)M(l,m) = M(l + l,m)L(l,m) 

and it is through this compatibility condition that we arrive at the integrable non- 
linear equation associated with the Lax pair. 

Despite the existence of a Lax pair often being used as the definition of integra- 
bility for a given equation [3], there have been few studies that sought to find or 
categorize nonlinear equations that used Lax pairs as their starting point. Of those 
studies that did begin with Lax pairs, most chose a form of the Lax pair a priori, 
that is an assumption was made concerning the dependence of the linear systems 
on the spectral parameter, thus limiting the possible results. 
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The present study begins with Lax pairs that are 2x2, where each entry of the 
Lax matrices contains only one separable term. The Lax pairs are otherwise general 
in that no assumptions are made as to the explicit dependence of any quantities 
within the Lax matrices on the lattice variables, I and m, or on the spectral variable 
n. By one separable term we mean that each entry contains a term that can be 
split into a product of two parts, one that depends on the lattice variables, and 
another that depends only on the spectral variable. For example, in the 11 entry of 
the L matrix, we write the separable term a(/, m)A(n). Both a and A may contain 
multiple terms themselves, say a = a%, A= Y) j Aj, however, all terms within a 
must multiply all those within A, the other entries their own similar term. So, we 
could have an L matrix 



where we have written the 11 entry in the expanded form and left the other entries 
abbreviated to save space, no other terms can be added into any entry. M must also 
contain just one separable term in each entry, although these terms are independent 
of those in L. 

The reason for limiting the Lax pairs to those that are 2x2 with one separable 
term in each entry is two fold. Firstly, Lax pairs with more terms typically lead 
to equations of higher order, as can be seen from hierarchies of equations with 
Lax pairs [U |2j , therefore we constrain our study to the lower order equations by 
limiting the number of terms. Secondly, we limit the number of terms present 
in the compatibility condition and thus render it less complicated to examine all 
of the combinations of terms that can arise there. A combination of terms in the 
compatibility condition defines a system of equations that we subsequently solve, in 
a manner that preserves its full generality, up to a point where a nonlinear evolution 
equation is apparent, or it has been shown that the system cannot be associated 
with a nonlinear equation. Testing all combinations of terms, we thereby survey 
the complete set of Lax pairs of the type described. 

In fact, of all the potential Lax pairs identified by this method, only two lead to 
interesting evolution equations. These are higher order varieties of the lattice sine- 
Gordon (LSG) and the LMKdV equations, which can be found in section ITTTl The 
remaining systems are shown to be trivial, overdetermined or underdetermined. As 
noted in [6] , this suggests a connection between the singularity confinement method 
and the existence of a Lax pair. 

As we do not make any assumptions about the explicit dependence on the spec- 
tral parameter, we show that a particular nonlinear equation may have many Lax 
pairs, all depending on the spectral parameter in different ways. The effect that 
this freedom has on the process of inverse scattering is, as yet, unclear. 

This paper is organized as follows: section ITTTl presents the major results, those 
being the higher order versions of LMKdV and LSG, as well as a statement of 
the completeness theorem. The method of identifying and analyzing the viable 
Lax pairs is laid out in section where a representative list of all the Lax pairs 
identified can be found. Section [3] explains how the higher order LSG and LMKdV 
equations are derived from the general form of their Lax pairs and sectionUJprovides 
examples that describe why most Lax pairs found in section[2]lead to trivial systems. 
A conclusion rounds out the paper. 

1.1. Results. Note that all difference equations in the remainder of this paper will 
utilize the notation 




x 



X 



x(l + 1, to), 

x(l, TO + 1). 
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As one of the two main results of the paper, we present two new integrable 
nonlinear partial difference equations. The first equation, 

LSG2: h \\H\xy — — — H — (1.2a) 

ax p x xy 



ay , A 2 /i 2 



XiPixy (l-2b) 



P 2/ xy cr y 

is referred to as LSG2 because it is second order in each of the lattice dimensions, 
and because setting x = y returns the familiar LSG equation, in a non-autonomous 
form. Here x — x(l,m) and y — y(l,m) are the dependent variables, Ai = Aj(£) and 

Pi = Pi(w) are parameter functions, as are p = A3 and a — p 3 . 
Similarly, the equation 

LMKdVV — — H — = A 2 cr- + ppi— (1.3a) 

a x p y y x 

ppixy + \ 2 axy — —xy-\ -xy (1.3b) 

p a 

where the terms are as for LSG2, is referred to as LMKdV2, again because setting 
x = y brings about LMKdV. Note that the LMKdV so attained is of a more general 
form than the most common variety listed in the beginning of the introduction. 
The second main result comes in the form of a theorem that we state as follows: 

Theorem 1. The system of equations that arise via the compatibility condition 
of any 2x2 Lax pair il.l)) with one nonzero, separable term in each entry of 
each matrix (as described in the introduction) is either trivial, under determined, 
over determined, or can be reduced to one of LSG2 or LMKdV2- 

The proof of theorem [T] lies in considering all of the possible sets of equations 
that can arise from the compatibility condition of such Lax pairs, and solving those 
sets of equations in a way that retains their full freedom. This proof occupies the 
remainder of the paper. 

The Lax pairs associated with LSG2 and LMKdV2 respectively are listed below, 
these Lax pairs are derived in section [31 The Lax pair for LSG2 is 

- - ( ^:,x 1% ) (-) 

M = ( F t l[GX) F ' /V ) (L4b) 
\ F 1 p 1 y F 2 a J 

While the Lax pair for LMKdV2 is 

l - ( F £f/: F rS: ) ^ 

M = ( Fl ^ lx F f /<T ) (1.5b) 
y F 2 ay F 1 p 2 ) 

where = Fi (n) are arbitrary functions of the spectral variable n with the condi- 
tion that F\ 7^ kF 2 , where k is a constant. 

2. Method of identifying potential Lax pairs 

From Lax pairs of the type we consider here, the compatibility condition produces 
a set of equations, each equation being due to one of the linearly independent 
spectral terms that arises in some entry. Studies that have searched for integrable 
systems by beginning with a Lax pair, whether an isospectral or isomonodromy Lax 
pair or otherwise, typically assume some dependence of the Lax pair on the spectral 
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parameter, then solve the compatibility condition for the evolution equation [H 
[TTJ HU] • Most often a polynomial or rational dependence on the spectral variable is 
used |13j . but any type of explicit dependence could be investigated, for example 
Weierstrass elliptic functions. 

Example 1. Consider the following L and M matrices 

L ( ap' Abp 

\ c dp' 

M = + & 

where p is the Weierstrass elliptic function in the spectral parameter n only, and 
all other quantities are functions of both the lattice variables I and m. From the 
compatibility condition, LM = AIL, and noting that p' 2 = jp 3 — §2p — 53 where 
gi = constant, we find the following equations in the 12 entry 



12 



p 3 : 


afl 


= ba-\ 


-dp 


P 2 : 


bS 


= 




P : 


a/3 


= dj3- 


ba 


P°: 


ap 


= dp 





(2.1) 



where we have separated out the equations coming from different orders of the spec- 
tral parameter. 

Clearly, this choice of L and M does not yield an interesting evolution equation 
through their compatibility, this example was instead chosen because it illustrates 
an important point. Multiplying together two functions of the spectral parameter 
can produce numerous orders that may or may not be proportional to other spectral 
term products in the compatibility condition. From example [TJ the spectral term 
multiplying p 2 did not 'match up' with other terms in the 12 entry at that order, 
which brought about the equation bd — 0, forcing some term to be zero. However, 
the other spectral terms did turn out more meaningful equations, highlighting the 
need to choose the dependence on the spectral parameter carefully so that none of 
the resulting equations force any lattice terms to be zero. 

The inclusion of a zero lattice term does not preclude the existence of an in- 
teresting evolution equation. However, we are essentially classifying Lax pairs by 
the number of terms in their entries and the class of Lax pairs presently under 
inspection contains one separable term in each entry of their 2x2 matrices. If any 
of those terms were forced to zero then the resulting Lax pair would actually come 
under a different category in the present framework. 

The most general form of a 2 x 2 Lax pair with exactly one separable term in 
each entry of the L and M matrices is: 

aA bB 



L = 
M = 



cC dD 



aA PE 

7 r 6A 

Where lower cases represent lattice terms and upper cases represent spectral terms. 
The compatibility condition is LAI = AIL, of which we initially concentrate on the 
12 entry. 

aPAE + bSBA = baBA + dpDE (2.2) 
At this stage we are only concerned with the various linearly independent spectral 
terms that appear. These will determine the set of equations that come out of the 
compatibility condition, which are subsequently solved to find the corresponding 
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evolution equation. That being the case, there are four quantities to contend with 
in this entry, AE, BA, BA and DE. Any of these four products can lead to multiple, 
linearly independent spectral terms, all of which must match up with at least one 
other spectral term from one of the three remaining products in this entry. If 
there exists some spectral term that does not match up with a spectral term from 
another product, then the lattice term that multiplies it will have to be zero, which 
is forbidden. 

Let us label the terms in each product as follows AE = J^i Fas ( > BA — J^i Fba ( > 
etc. All the spectral terms that occur in the 12 entry of the compatibility condition 
can be sorted into four groups according to the lattice terms that they multiply. 



12 



(a/3) (bS) | (ba) (dp) 

Fa-! FbAi FbAi Fnsi 

Fa~ 2 F BA2 F BA2 F DS2 ( 2 - 3 ) 



Where the line separating the four groups marks the position of the equals sign in 
the associated lattice term equations. Organizing the terms from example [T] in this 
way leads to 



12 



(a/3) (b5) (ba) (dp) 

P 3 P 2 P 3 P 3 

P I P P 

P° I P° 

The present study will utilize the word 'group' in the general English sense, our 
groups refer to collections of spectral terms that multiply the same lattice terms in 
an entry of the compatibility condition. We have already seen that all terms in all 
groups must be proportional to a term from at least one of the other three groups 
in this entry. Conversely, where there are spectral terms that are proportional to 
others from other groups, an equation relating the corresponding lattice terms will 
thus be defined. Still with the above example, that Lax pair has a dependence on 
the spectral parameter such that the groups multiplying a/3, ba and dp all contain 
the spectral term p 3 , for which the corresponding equation in the lattice terms is 
a/3 = ba + d/3. 

At this stage, the number of possible sets of proportional spectral terms, and 
therefore the number of possible combinations of equations yielded by the compati- 
bility condition, is unmanageably large. We require further considerations to bring 
the problem under control. 

2.1. Links and equivalent equations at different orders. The following defi- 
nition applies to sets of spectral terms within an entry of the compatibility condition 
that bring about the equations used to find the evolution equation. The word pro- 
portional is used in the sense that two terms F\ and F^ are proportional if F± = kF^ 
for some finite constant k. 

Definition 1. A link is a set of proportional spectral terms in the same matrix 
entry of the compatibility condition. A set of two proportional terms is a single 
link, three terms a double link, and four terms a triple link. 

Naturally, the spectral terms that comprise a link must each reside in a different 
group of spectral terms within an entry. If there are two or more spectral terms 
that are proportional to one another within the same group, they are simply added 
together to make one term. Since each group of spectral terms multiplies the same 
lattice term, one may speak of either links between the groups of spectral terms or 
links between lattice terms, with the same meaning. The above definition captures 
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the idea that the entries of the compatibility condition give rise to different lattice 
term equations at different orders in the spectral term, without appealing powers 
of some basic function. 

With the employment of proportional terms comes the possibility of constants 
of proportionality, which we begin to deal with here. 

Fact 1. If there exist two distinct single links between the same two groups of 
spectral terms, then the corresponding constants of proportionality must be equal. 

The proof of this fact is elementary: say that one single link is formed by allowing 
FaSx oc -Fbsi , where neither of these terms is proportional to any other spectral 
term that arises in the 12 entry, and the other single link corresponds to Fas 2 (x 
Fde.21 where again these terms link with no others. By including some constants 
of proportionality, k\ and ki respectively, we can write down the equations that 
correspond to these two single links, those being a/3 = k\d(3 and d[3 — k2d(3. Since 
both equations must hold, and none of the lattice terms can be zero, it is clear that 
we must have k\ — k%- 

This simple fact proves to be rather important because it ensures that all links 
between the same two groups can be bundled together. Further, all the spectral 
terms, in some group, that correspond to those single links with one other group, 
can be treated as a single spectral term. So, where Fas, 1 and formed one single 
link and Fab. 2 an d Fos 2 formed another between the same two groups, we can lump 
together Fasi + Fas 2 = Gi an d be sure that it links with Fob.! + Fds 2 — kG\. 

Still with the 12 entry, if there exist multiple double links between the same 
three lattice terms, a(3, ba and d[3 say, then the lattice term equations that result 
from those links can be written 



Where K is a matrix of the constants of proportionality between the various spectral 
terms, normalized so that each entry of the first column of K is unity. Each 
row of K corresponds to a double link. If K is such that equation (|2.4|) is over- 
determined or uniquely solvable, then any Lax pair possessing the corresponding 
links is inconsistent or contains a zero lattice term. Therefore, we need not consider 
more than two double links between the same three spectral terms, although there 
may be multiplicity within those two double links. By the same argument we can 
allow a maximum of three different triple links between the same four lattice terms 
in an entry of the compatibility condition. 

2.2. Link symbolism. The abundance of Lax pairs that need to be checked ne- 
cessitates the introduction of a shorthand, which will be based on their links. The 
off-diagonal entries both contain four groups of spectral terms, each of which mul- 
tiplies a single product of lattice terms. The 12 entry contains the spectral term 
products AE, BA, BA, and DE associated with the lattice term products af3, bS, 
ba and d[3 respectively. For the shorthand, we always set out the spectral terms in 
the same way on the page 



Each link can be represented by lines between the quantities that are proportional 
to each other, and so we will use the symbols listed in table [1] to represent the 
combinations of links in the 12 entry, where Fas is some term from the group of 
spectral terms formed by taking the product AE, and other terms are similarly 
labeled. 




(2.4) 



AE BA 
BA DE 
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SvTnhnl 


Links 


X 


Fas oc F ds , F B a oc F ba 




Fas oc Fba, Fba oc Fde 


□ 


Fas oc F BA , Fba oc F dh 


Fas oc F ds cx F ba oc F B a 


171 


Fas oc F ba oc F BAi , Fba^ oc F ds 



Table 1. Symbols used to represent the link combinations in the 
off-diagonal entries. Note that Fba* ^ kFsAj , k a constant 



The 21 entry is similar to the 12 entry in that it possesses four distinct products 
of spectral terms. The same symbols listed in table Q] are used again for the 21 
entry, with clear meaning given that the spectral products are set out as follows 

CA AT 
DT CA 

However, the diagonal entries are slightly different because each diagonal entry 
contains a product that occurs twice: AA occurs twice in the 11 entry and DA 
twice in 22. This automatically causes the associated lattice terms to be paired in 
their respective entries and, as such, AA and DA need not be linked with another 
spectral term to prevent a zero lattice term. This being the case, there are really 
only three spectral term products to consider in both of the diagonal entries, one of 
which need not be linked to the other two, and our symbols reflect that. Positioning 
the spectral term products as follows 



AA 



BY 
CE 



We symbolize the links as indicated in table [21 where the symbol '< — ' is used to 
represent the repeated spectral term products. 



Symbol 


Links 


~< 


Faa oc F B t oc F C s 


~> 


Faa alone, F B t oc F C s 


1 ' 


FaA; oc Fbt and Faa^ oc Fcs separately 



Table 2. Symbols used to represent the link combinations in the 
diagonal entries of the compatibility condition. Note that Faa 4 ^ 
/cFaAj j k a constant 



2.3. Which Lax pairs need to be checked? All link combinations that do not 
force any lattice terms to be zero are checked systematically. The procedure for 
doing this runs as follows: 

(1) Begin with the 12 entry, assume that only single links exist there, and list 
all single link combinations that produce different lattice term equations. 
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(2) Construct the links found in the previous step by choosing proportional 
sets of spectral terms in the appropriate groups. 

(3) Move to the 21 entry, noting the spectral term constructions from the pre- 
vious step, and identify all the viable link combinations in this entry. 

(4) Repeat the previous step in the diagonal entries. 

After identifying all link combinations with single links in the 12 entry, we repeat 
the entire process assuming that double and possible single links exist there, and 
then repeat once more with triple, and possibly double and single links in the 12 
entry. Finally, the corresponding set of lattice term equations for each Lax pair 
must be analyzed to find the resulting evolution equation. This analysis needs to 
be conducted in a manner that preserves the full freedom of the system, as described 
in section [3] 

2.3.1. Single links in the 12 entry. The four groups of spectral terms in the 12 entry 
are 

AE BA 
BA DE 

Each group must be linked to another so, using single links, the group AE must be 
linked to at least one of the three other groups, BA, BA, or DE. For arguments 
sake, say that there exists a single link between AE and BA. 

AE BA 
BA DE 

That leaves both groups BA and DE requiring links and, since we are only con- 
cerned with single links at the moment, these two groups of spectral terms can be 
linked to each other, or to one of AE or BA, in distinct, single links. However, it 
can be shown that if any group in cither off-diagonal entry possesses single links 
between it and two other groups, the resulting Lax pair is associated with a trivial 
evolution equation (see proposition Q] below) . Hence, given the first link between 
AE and BA, the only other single link that needs to be considered is between BA 
and DE. Therefore, table [3] lists the only single link combinations in the 12 entry 
that require further analysis. 



AE BA 
BA DE 



12 



Table 3. Single link combinations in the 12 entry 



Proposition 1. // there exist two single links between some lattice term and two 
others in an off- diagonal entry, then the resulting evolution equation is trivial. 

The proof of Proposition [1] lies simply in checking all the possible link combina- 
tions that meet the criterion. In that sense there is nothing special about Lax pairs 
of this type, they are only treated separately because there are many such cases, 
so including them with the others would unnecessarily lengthen the argument. 

Before links in the 21 entry are examined, we must select spectral terms that 
engender the links already chosen in the 12 entry. We shall proceed with the 
analysis under the assumption that single links between AE and -BA, and between 
BA and DE, although the other combinations can be dealt with in the same way. 
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The specified links are constructed by choosing spectral terms 

A = -(Fi + ...) (2.5a) 

D = 1(F 2 + ...) (2.5b) 

A = + (2.5c) 

A = i(F 2 + ...) (2.5d) 

where F{ — Fi(n) and Fi ^ kF 2 , k a constant. It is understood that while there is 
room for other spectral terms in the expression A = =(Fi + ...), there cannot be a 
term proportional to F 2 /'E. in A, as this would cause the single link corresponding 
to F2 to become a double link. In this way the desired links are constructed, 
plus we have allowed for additional links should they be appropriate or required 
later. Readers may note the omission of any constants of proportionality in the 
above, however, the constants that could have been written at this point can all be 
absorbed into the lattice terms that they multiply. 

Turning our attention to the 21 entry of the compatibility condition, caCA + 
d"fDT = 07 AT + cSCA, the following spectral terms appear 



21 



(ca) (dj) I (07) (cS) 
c 771 r tp 1 r 7? c tp 



5^2 (2.6) 



Since every spectral term must be proportional to another in the same entry, 
jjFi must be proportional to =i*i or =F 2 . Clearly, %F\ cannot be proportional 
to §1*2) nor can it link with some other term that we are yet to define, as this 
would introduce single links between some lattice term and two others, the situation 
excluded by proposition [TJ Moreover, we can exclude the case with §Fi x ^F 2 
because this also requires the remaining spectral terms to be proportional to one 
another, i.e. %F 2 oc =i*i. These two conditions on the spectral terms imply that 
F\ k F2, contradicting a previous assumption. Hence, the links chosen in the 
12 entry leave only one choice for the links in the 21 entry of the compatibility 
condition: j^Fx x =Fi and ^F 2 x =F 2 , which can be written more succinctly as 
CZ = BT 



12 AS BA 21 CA AT 



BA DE DT CA 

The terms arising in the diagonal entries are listed below 

m\ 1 

(da -ad) F?/(BZ) \ (c/3 - 67) BT 

(2.7) 

I 

(dS-dS) Fi/(BE) I Qyy-c0) BT 

Note that the terms in equation (|2.7|) are set out slightly differently to those in 
table because we have already determined that BT x CS. 

The spectral terms in the diagonal entries, in this case, do not necessarily have 
to link with others because they are multiplied by more than one lattice term, as 
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indicated in equation (|2.7p where the lattice terms appear in parentheses to the left 
of the spectral terms they multiply. A lone spectral term in the diagonal entries, 
given the links already constructed in the off-diagonal entries, will not bring about 
a zero lattice term. Consequently, there exist three possible link combinations for 
the diagonal entries: BY oc F^/(BE), BY oc F$/(BE) or BY is proportional to 
neither F 2 / (BE.) nor F± / (BE). All three choices lead to trivial evolution equations 
and we shall continue the analysis under the assumption that BY oc F^/(B3). A 
link combination has now been chosen in each of the entries of the compatibility 
condition, these are shown in table 0] 



12 


21 


\m 


22 


AS BA 


CA DT 


BT 

AA 


BF 

DA 


BA DS 


AT CA 


CS 


CS 








l | 

' 1 \ 



Table 4. An example of the links that define a Lax pair 



Gauge transformations can be used to remove the dependence on some of the 
spectral terms and, as such, we expect some redundancy. The links constructed 
above are achieved by setting the values of the spectral terms to 



A 


= F 1} 


A 


= F 1 


B 


= 1, 




= 1 


C 


= Fl 


r 


= F? 


D 


= F 2 , 


A 


= F 2 



where F\ and F 2 are any functions of the spectral variable n, such that F\ ^ kF 2 , 
k = constant. Note that one of F\ or F 2 may itself be a constant. No constants 
of proportionality are required as these can be absorbed into lattice terms in this 
case. The same links are reproduced by any suite of spectral terms that satisfies 
the following conditions, brought about by the links described above 

A = Fx/3, A = Fi/B 

C = if/E, T = Ff/B 

D = F 2 /E, A = F 2 /B 

The resulting set of equations that are produced by the compatibility condition 
for this Lax pair are written in (|2.8p . although they can be read from the links 
given in tabled] 



aa 


— aa = 


c/3 — &7 


dS 


-dS = 





b", 


-c/3 = 







a/3 = 


ba 




u = 


dp 




ca = 


aj 




dj = 


cS 



This rounds the description of the Lax pairs with only single links in the off- 
diagonal entries, in practice one would continue by analyzing equations (|2.8j) to find 
the associated evolution equation, which in this case is trivial. On the possibility of 
including extra spectral terms to augment the links used here, see equation (|2.5p . 
we note that fewer terms allow greater freedom and that any additional links could 
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only lead to a more constrained system, one that certainly could not sustain an 
interesting evolution equation considering that the less constrained example here 
leads to a trivial result. Also, the alternative links that cause there to be one 
equation in the 22 entry and two equations in the 11 entry of the compatibility 
condition, see after equation (|2.7p , lead to the same evolution equations found here 
as the Lax pair is symmetric in that sense. 

2.3.2. Double links in the 12 entry. Here link combinations that consist of double 
and possibly single links in the 12 entry are investigated. There are more possi- 
bilities in this class than when only considering single links, however the number 
is reduced by noting that some sets of link combinations are equivalent from the 
perspective of the lattice term equations. For example the following pair of link 
combinations in the 12 entry are clearly equivalent. 

ae ba rz ^ 

BA DE — 



12 



Also, all link combinations that possess two double links in the 12 entry are nearly 
equivalent, enough to consider them all together. The equivalence is because the 
lattice term equations corresponding to any two double links in this entry can be 
manipulated so that the they are the same as those from any other combination of 
two double links. The difference that may arise comes from the diagonal entries, 
where the particular pair of double links chosen in the 12 entry can affect the variety 
of link combinations possible. However, the difference is not sufficient to alter the 
overall outcome that these systems are overdetermined. The complete list of double 
link combinations in the 12 entry is listed in table [5] 

A3 BA rZ ^ V \ N ^ ^ n 
BA DE 



12 



Table 5. All double link combinations in the 12 entry that require analysis 



To exemplify the method of construction of Lax pairs with double links in the 12 
entry, we choose link combination in table [5] where there is a double link between 
AE, BA and BA and a single link between BA and DE. Spectral terms that 
generate this choice of links are given in equation 



A = F 1 /E, A = Fx/B 

A = (Fi + kF 2 )/B, D = F 2 /E [Z ^> 

where our usual nomenclature applies, i.e. Fi — Fi(n), k = constant. 

Using the expressions found in the 12 entry, the terms that arise in the 21 entry 
are F\C/B, F 2 T/E, FiT/E and F 2 C/B, where FiC/B arises twice. It is convenient 
to rearrange the terms found in the 21 entry into columns of terms that cannot be 
linked, this is done in equation (|2.10[) . 

F 1 C/B FiT/E F\T/E F lC /B~ 

F 2 T/E {F 1 + F 2 )C/B F 2 T/E F 2 C/B ^ iUj 

F 2 T/E must link with F 2 C/B because linking with FiC/B would lead to the 
contradictory F-y oc F 2 , since that would also necessitate a link between F 2 C/B 
and FiT/E. That leaves two possibilities: BY oc CE, or we can split F 2 C/B oc 
+ F 2 )T/E, noting that the multiplicity of the term F\C jB means that it need 
not link to another in this entry. The second possibility is neglected, though, as 
it gives rise to a zero lattice term in one of the diagonal entries to be considered 
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below. Thus, by a process of elimination, the links selected in the 12 entry leave 
only one choice for the 21 entry, which is shown in equation (|2.1ip 



12 AE BA 

V 

BA DE 



21 CA 
DT 



AT 
CA 



(2.11) 



Move now to the 11 entry where the relevant spectral term products are AA oc 
Fi I (BE), and BY oc CE, see equation (|2.12p . We choose to link these two spectral 
terms to get one equation in the 11 entry and find that there are necessarily two 
equations in the 22 entry. Note that the opposite case where there are two equations 
in the 11 entry and one in the 22 entry can also exist by splitting up BY into two 
terms, however this will lead to the same evolution equation. 



m 



[act - aa) F?/(BE) 



22 



(c/3 - 67) BY 



(2.12) 



F X F 2 /(BE) 

(dS-d5) I (67 -c/3) BY 

Fi/(BE) I 



Using our symbolism, the links that define this Lax pair are shown in table El 



12 


21 


m 


22 


AS BA 


CA DT 


BT 

AA 


BF 

DA 


BA Dr. 


AY CA 


CE 


CE 


V 




-< 





Table 6. Links for a Lax pair with double links in the off-diagonal entries 



The spectral term relations that must be satisfied to bring about the links in 
tableware given in (|2.13[) . 



A = 


Fi, 


A 


Fi 


B = 


1, 




= 1 


C = 


Fl 


r 


= F? 


D = 


F 2 , 


A 


= Fx + l 



(2.13) 



where F\ and Fi are any functions of the spectral variable n, such that F\ ^ kF 2 , 
k = constant. Notice that a constant appears in the expression for A in equation 
(|2.13[) . since one of the constants of proportionality in this expression cannot be 
absorbed into the multiplying lattice term, 5. 
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The lattice term equations that arise via the compatibility condition from this 
Lax pair arc shown in equation (|2.14p . 



da — 


aa — 


c/3- 


\r t 


dS - 


dS = 







67 — 


c(3 = 







a/3 4 


-bS = 


ba 






b5 = 


kdp 






ca = 


c?7 + 


cS 




07 = 


kcS 





We thus conclude the description of the formation of Lax pairs with at least one 
double link in the 12 entry. Naturally there are other Lax pairs of this type but 
they are formed in a similar manner to that described here. 



2.3.3. A triple link in the 12 entry. Lastly, link combinations including triple, and 
possibly double and single links in the 12 entry must be considered. It is not difficult 
to see that the only possibility that needs to be investigated is that with one triple 
link between all four lattice terms in the 12 entry, any additional links in this entry 
constrain the problem too heavily and lead to trivial evolution equations only. 

A triple link in the 12 entry of the compatibility condition can be formed by 
setting the spectral terms to those in equation (I2.15|) below. 



A = F 1 /E, A = Fx/B 

A = F 1 /B 1 D = F 1 /E 



12 



AE BA 

□ (2.15) 
BA DE 



Given the values in (|2.15p . the resulting spectral terms that appear in the 21 
entry are as shown in equation (|2.16[) . 



21 



CA AT FrC/B F 1 Y/E 

DT CA F 1 T/E F 1 C/B l ^ iDJ 

We therefore have two possibilities in the 21 entry depending on weather or not 
BY oc CE. If BY is proportional CE then there is another triple link in the 21 
entry, otherwise there is a pair of single links. The latter case yields the links given 
in item 4 of tableland leads to a trivial evolution equation, while the former yields 
Lax pairs that include one for the LMKdV2 system. 

We continue the analysis here assuming BY oc CE, in this case the spectral terms 
in the diagonal entries are shown in equation (|2.17[) . 



22 



{aa-aa) F%/{BE) | (c/3 - 67) BY 

(2A7) 

I 

(dS-dS) F?/(BE) I (67-C/3) BY 



Again, two possibilities present themselves, this time depending on whether F? /(BE) oc 
BY. The case where the proportionality does not hold provides a Lax pair for the 
LMKdV2 system, which is discussed in section [3] When F^ /(BE) oc BY does hold, 
a unique situation unfolds where each entry of the compatibility condition contains 
only one equation. This special case is discussed in section [ 
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2.4. List of link combinations. We are now in a position to tabulate the results. 
Table [7J contains a representative selection of all possible link combinations for 2x2 
Lax pairs with a single, separable term in each entry of the L and M matrices. 
There are still other link combinations that were analyzed but do not appear in 
table [7] because they are equivalent to a combination that does appear, or because 
it is clear that the corresponding Lax pair cannot yield an interesting evolution 
equation since a similar, less constrained combination of links is listed as trivial or 
over-determined. 





1 9 
1Z 


91 
Zl 


rm 


99 

zz 






AE BA 


CA DT 


bf 

AA 


Br 

DA 


Evolution Eqn 




BA DB 


AT CA 
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X 


LSG 2 (IL2j) 
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n 
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LMKdV 2 (|1.3p 


3 


v 
/\ 


y 


/ 


/ 


Trivial 


4 


V 

/x 


n 


1 

>— 1\ 


i-is 


Trivial 


5 


y 

/X 




1 | ^ 

1— 1\ 


>— '•v 


Zero term 


6 






m/ 

X 


/ 


Trivial 


7 






^/ 




Underdetermined 


8 


V\ 


N 




~< 


Trivial 


9 


V 




-< 




Trivial 


10 


\K 




~^ 


-< 


Overdetermincd 


11 


N 


171 




-< 


Trivial 


12 








-< 


Overdetermined 


13 


A 






-< 


Overdetermincd 


14 


n 






-< 


Overdetermined 


15 


□ 


□ 


-< 


-< 


Special case 



Table 7. List of link combinations used to construct possible Lax pairs 



Items 1 and 2 in table [7] are analyzed thoroughly in section [3J The other represen- 
tative link combinations are dealt with in section fU 

3. Derivation of higher order LSG and LMKdV systems 

While fourteen potentially viable types of Lax pairs were identified in section 
[3J only two types lead to non-trivial, well determined evolution equations. The 
two systems thus found are LMKdV2 and LSG2 and this section describes the 
derivation of these systems from the general form of their Lax pairs. It is important 
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to note that no freedom in the lattice terms is lost through this process, all values 
that the lattice terms take are dictated by the sets of equations effectuated by 
the compatibility condition. As such we conclude that the systems so derived (or 
equivalent systems) are the most general ones that can be associated with their Lax 
pairs. 

In fact, such calculations have been performed many times before and yet neither 
LSG2 or LMKdV2 appear to have been published, despite coming from Lax pairs 
with simple forms that have certainly already been considered previously H2] . 
Hence, it is necessary to outline the method used to derive LSG2 and LMKdV2 in 
detail. 

3.1. LSG2. For LSG2 the Lax pair used as a starting point is 

* - ( % S ) 

" - ( % Is ) 

Where a, b, c, d, a, (3, 7 and 6 are all functions of both the lattice variables I and 
m, referred to as lattice terms, and F\ and F? depend on the spectral variable only, 
we refer to these terms as spectral terms. 

It is possible to remove some of the lattice terms using a gauge transformation. 
While this would reduce the complexity of the system, it is not clear at this point 
which of the lattice terms would best be removed. Experience shows that natu- 
ral transformations present themselves in the course of solving the compatibility 
condition and so we shall wait until later to remove some lattice terms, keeping in 
mind that we expect some freedom to disappear from each of the Lax matrices L 
and M. 

The compatibility condition for (|3.ip . LM = ML, leads to the following system 
of difference equations in the lattice terms 

(3.2a) 
(3.2b) 
(3.2c) 
(3. 2d) 
(3.2e) 
(3.2f) 

Some of equations (|3.2p are linear and some nonlinear. The linear equations can 
be solved easily when in the form 

ki4> — ^20 = k^ip — kiifi (3-3) 
where <f> and ip are lattice terms and fcj constants. Equation (|3.3[) implies 

if) = k\v - k 2 v + fj,(— y 

- k 3 v - k iV + A(^) m 
ki 

where we have introduced the new lattice term v — v(l,m), and A = A(Z) and 
/i = /Lt(m) are constants of integration. The same fact also applies in a multiplicative 
sense, in particular 

d> ib v v 



aa + 67 = 


aa + cf3 


dd + c/3 = 


dS + b^ 


a/3 = 


dp 


is = 


ba 


ca — 


cS 


c?7 = 


07 
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We now proceed to solve the system (13.21) . Equations (|3.2cj) to (|3.2ij) are linear 
and may solved in pairs as follows. Multiply equations (|3.2cp by (|3.2fj) to find 
ad/ (ad) = P^/(P^), which is integrated for 

a = \lv 5 /(v 5 d) 

P = M5W 5 /(«57) 

where we have introduced Ao = ^o(l), P5 = Ms( m ) an d V5 = vs(l, rri). Now use 
equation (|3.2c|) again to find v^/{v^) = j-j^ and integrate for 

d = \opv 2 /v2 

where p — ^(l)^ 1 * 1 . Substitute these values back into the expressions for a and 
f3, and replace v$/v2 1— ► V\, and M5/M2 1— » Mi resulting in 



(3.4a) 



d = A p— (3.4b) 

13 = vi— (3.4c) 

7 = A*2— (3.4d) 
Perform similar calculations on equations (|3.2d[) and (|3.2ep to find 

b = Ai— (3.5a) 

Vi 

c = A 2 — (3.5b) 

a — — 

a v 3 

S = Pqct— (3.5d) 



(3.5c) 



Where a = ^(m)^ -1 ^ . 

When equations (13. 4|) and ()3.5|) are substituted into (|3.2a|) and (|3.2bj) we find 
the following equations respectively 

, MO V!V 3 V 2 V 3 A ,„ „ v 

AoP — h A1M1 — — = — Mo<7 — — + A2M2- (3.6a) 

— Moc- V \2P2- = A p — + A1//1- — (3.6b) 

p V2V4 V3V2 o V2V4 V1V4 

Multiplying (|3.6ap by Ui/©i and (|3.6b|) by V4/V4 indicates that certain variables 
always appear in combination. As such, we set Vx, V2 = 1 without loss of generality, 
and rename V3 = x and V4 = y. This is the manifest reduction in freedom that was 
expected from the perspective of gauge transformations. The parameter functions 
similarly appear in ratios, hence we set Ao = 1 and po = 1 without loss of generality. 
Making the substitutions we arrive at a pair of nonlinear partial difference equations 
in x and y, with arbitrary non-autonomous terms Aj(Z) and pi(m), which together 
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form LSG2, the evolution equation associated with the Lax pair. 

px . ox X 2 p 2 ,„ _ x 

h Mpixy = — — H — (3.7a) 

ex x p x xy 

<J§ . A 2 /i2 PS . . ,„ „ s 

— -H ; — = hAi/xixy (3.7b) 

py xy ay 

This pair of equations can be thought of as a higher order lattice sine-Gordon 
system because the lattice sine-Gordon equation (LSG) is retrieved by setting y = x 
in either expression 

LSG : xx( Ai/iixx) = — xx — X2P2 

P v 

The Lax pair Q1.40 for LSG2 is obtained by substituting the calculated values of 
the lattice terms back into (13.11). 



3.2. LMKdV"2. The the general form of the Lax pair for LMKdV2 is similar to 
that for LSG2, the only difference being that here both matrices exhibit the same 
dependence on the spectral variable, while the dependence was antisymmetric in 
the previous case. 

' - ( Z Z ) <") 

»-{% % ) <-> 

As above, a, b, c, d, a, 0, 7 and 8 are all functions of the both the lattice variables 
I and m, F± and F% depend on the spectral variable only. 

The compatibility condition leads to six equations coming from the different 
orders of the spectral variable in each entry. 



aa = 


a a 




(3.10a) 


67 = 


cfi 




(3.10b) 


dS = 


dS 




(3.10c) 


c/3 = 


&7 




(3.10d) 


d(3 + b5 = 


ba ^ 


-dp 


(3.10e) 


COL + dj = 


a7 - 


- c5 


(3.10f) 



Equations p.lOap and (|3.10cp are integrated immediately, while (|3. 10b|) and (|3. 10d[) 
are multiplied together and dealt with by a similar method to that used for the LSG2 
system of section 13.11 leading to the following results 



a = 


Xivi/vi 


(3.11a) 


b = 


X pv 3 /vi 


(3.11b) 


c = 


Xov 4 /(pv 3 ) 


(3.11c) 


d = 


X2V2/V2 


(3. lid) 


a = 


Hivi/vi 


(3. lie) 


P = 


p. V3/(av 4 ) 


(3.11f) 


7 = 


p,oov 4 /v 3 


(3-llg) 


S = 


V2V2/V2 


(3.11h) 



Where Vi = Vi(l,m), Xi — Xi(l), pi = pi(m), p ~ X 3 ^ and a = p 3 ^ . Substitut- 
ing these values into the two remaining equations, (|3.10e[) and (|3.10fj) . shows that 
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terms consistently appear in ratios, as they did with the LSG2 system. We choose 
to set V2 = t>3 = 1, Ao = /Xo = 1) V\ = x and V4 = y and arrive at LMKdV2: 

\\x a 2 y . y . x 

— — H = A 2 cr-+p^i- (-J-12) 

a x p y y x 

ppixy + \20xy = —xy-\ -xy (3.13) 

pa 



4. HOW MOST LINK COMBINATIONS LEAD TO BAD EVOLUTION EQUATIONS 

Here we explain how most of the potentially viable Lax pairs fail to produce 
interesting evolution equations. The topic is split into four parts dealing with Lax 
pairs that lead to trivial, over-determined, under-determined evolution equations 
and the special case of item 15 in table [7J 



act 


= aa, 


67 = 


c3 


dS 


= dS, 


c/3 = 


717 


a(3 


= W, 


bS = 


ba 


07 


= ai, 


ca = 


cS 



4.1. Lax pairs that yield only trivial evolution equations. This section is 
pertinent to items 3, 4, 6, 8, 9 and 11 in tabled those Lax pairs who's compatibility 
conditions can be solved to a point where only linear equations remain, or the 
equations can be reduced to first order equation in one lattice direction only. 

The simplest route to triviality is to have a Lax pair with a set of equations that 
are all linear, as per item 3 in table [7] The compatibility condition leads to a set 
of eight linear equations in the eight initial lattice terms. 



(4.1) 



Equations ()4.ip can be solved easily using techniques described in section^ with 
the result of a simple, linear evolution equation. However, it is not necessary to 
conduct such analysis on this system since all equations (|4.ip are linear and they 
can not be expected to produce a nonlinear evolution equation. For this reason, 
other examples of link combinations that produce only linear equations have been 
omitted from table [7J 

Item number 8 from table [7J is an example of a Lax pair that leads to a trivial 
evolution equation in a more complex way. The equations that come out of its 
compatibility condition are 

(4.2a) 
(4.2b) 
(4.2c) 
(4.2d) 
(4.2e) 
(4.2f) 
(4.2g) 

Where k is a constant. This system of equations is solved as follows: integrate 
(I4.2a[) for a = Xv/v, a = fJ-v/v, introducing A = A(^), fi = p(m) and v = v(l,m). 
Multiply equations (|4.2bp and Q4.2g[ ) , then divide the product by (|4.2ep and use the 
values calculated for a and a to find 





a a 


= aa 






= c/3 


dS- 


f c/3 


= dS + b^ 


a/3 


+ bS 


= ba 




kba 


= -dp 




ca 


= aj + cS 




kca 


= — c?7 
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Equation (|4.3[) is integrated for b, which is substituted into equations (|4.2ep and 
(|4.2bj) yielding the following values 





dv 


b = 


\ 2 — 




cv 






= 






cv 




cv 


7 = 


— K/i — 




dv 



There are now three equations left to solve, (|4.2c[) . (|4.2d[) and (|4.2fj) . where it is 
found that c and v always appear as a product which suggests we introduce u = cv. 
Equation (j4.2c|) is then used to write 

6 = fJL h K/.1A — 7 

u ud 

Finally, we introduce x = ud/u and achieve the final two equations in x 

\ 2 (x-x) = fc(A 2 A-AA 2 ) (4.4a) 

x-x + k(\-\ 2 ) = k^(n\-\ 2 ) (4.4b) 

x 

Equations (|4.4|1 is an overdetermined set of two equations in the one variable, x, 
which cannot hope to yield an interesting evolution equation. This is because (|4.4aj) 
can be used to remove the dependence of x on m (m is the independent variable of 
the ' " ' direction), leaving x with a first order dependence on I in (|4.4bj) at best. 

4.2. Overdetermined Lax pairs. Here items 10 and 12 through 14 in table[7]are 
dealt with. Such Lax pairs have compatibility conditions that boil down to more 
equations than there are free lattice terms. These are not necessarily trivial, as 
the solution to one equation may solve another as well, that is it may be possible 
to make one or more equations redundant. Or one equation may be a compati- 
ble similarity condition for another equation as explained in [23j . Some systems 
of this type, where any hope of supporting an interesting evolution equation has 
been quashed, have already been considered in section 14.11 The remaining over- 
determined systems are considered here, however we do not attempt to resolve the 
issue of whether these systems support interesting evolution equations, we simply 
list them as being overdetermined. 

The most interesting instances of this type of system arise from Lax pairs with 
two double links in the off-diagonal entries, those represented by item 14 from table 
[7] All such systems, with minimum constraint on the lattice terms in the diagonal 
entries of the compatibility condition, lead to the same evolution equations. An 
example of a Lax pair with two double links in the off-diagonal entries is 

Fia b 
(F? + kiFiF 2 )c F 2 d 

M { (F? + feJifah (Fx + hF2)6 ) (4 ' 5bj 

where lower case letters except fcj are lattice terms, ki are constants of proportion- 
ality and Fi are spectral terms with F± not proportional to F 2 . This particular 
Lax pair is especially interesting because the solution of its compatibility condition 
involves integrating both additive and multiplicative linear difference equations, as 
described near equation ()3.3p in section [3] The final evolution equations achieved 
are relatively complicated and nonlinear, however there are two equations for one 
variable and it remains to be seen whether they can be reconciled. The following 
outlines how to solve the compatibility condition. 



(4.5a) 
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In the compatibility condition, LM — ML, the 11 entry dictates ki — ky, = k^. 
Redefine Fi i— > kiFi to remove k 2 from everywhere except the 22 entry of L where 
we set fco = 1/ K% and thus arrive at the following equations from the various orders 
of the compatibility condition 



aa + frj 


— aa + c/3 


(4.6a) 


c/3 


— fry 


(4.6b) 


d6 


= dS 


(4.6c) 


a/3 + bS 


= ba 


(4.6d) 


kibS 


= ba + kod(3 


(4.6e) 


ca 


= aj + cS 


(4.6f) 


ca + kgdj 


— k\c6 


(4.6g) 



Integrate (|4.6cj) for d — Xv/v, S — fiv/v, introducing v = v(l,m), A = X(l) and 

fi = jU(m), and use equation (|4.6b[) to see that (3 — fry/c. Now define s = c/v, 
t = r y/v and u — bv so that 

(3 = bi/s (4.7a) 

P~6gl) a=~(k 1 us-k Xt) (4.7b) 
s 

([4~6f)l a=l((k 1 -l)fis-k Xt) (4.7c) 

Substituting these (|4.7a|) values into equation (|4.6ep and rearranging leads to 



^) Mi) Hi) Mi) (4 - 8 ' 

which is an additive linear difference equation that can be integrated to find 

us 

= w~w + X 2 (4.9a) 

XX 

ui . . 

— = w — w + fj,2 (4.9b) 
Xfi 

The two equations that remain, equations (|4.6a|l and (|4.6d|) respectively, are 
written in terms of s, t and u 

s t s ~ s — t =s 

■±((k — l)ji — ko\~)(kiu fc A) + us = {(k — l)u fc A)(fci^ - fc A-) + ut- (4.10a) 

I s t t S t 

s t us s — 

-((fc-l)/t-fc A T )+^-= = fei/i^-fcoA (4.10b) 
I s ut t 

Make a further change of variables x = ui, y = us, to completely remove u from 
equations (|4.10[) and (|4.9p . Equation (|4.10b|) becomes 

//| + - 1)A - k a X*) = k^ - fc A (4.11) 

xx y x 

while equation (|4.10bj) becomes 

?((Al-l)A-*dAT)(*lM--M)+» = ((fci-l)A*--*»)A)(feiM-M*)+S- (4-12) 
x y x x y x 

where x and y are both written in terms of w according to 

x = — —(w — w + (12) (4.13a) 

y = —(w-w-\ 2 ) (4.13b) 
fci 
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Hence, there are two complicated, nonlinear equations for the one variable w, (|4.11[) 
and (14.12j) . These two equations may or may not be reconcilable, one way that 
they could be reconciled is if one equation was shown to be a compatible similarity 
constraint for the other [22l [23] 

4.3. Underdetermined Lax pairs. There are Lax pairs who's compatibility con- 
dition can be solved completely, while still leaving at least one lattice term free. In 
these cases there is no genuine evolution equation, although the freedom inherent in 
the system could to cause it to appear as though there was. In fact, any evolution 
equation, trivial, integrable or even chaotic, could be falsely represented by such a 
Lax pair. One such case is item 7 on table [7J that has as its compatibility condition 



da + by = 


aa + cj3 


(4.14a) 


dS + c(3 = 


dd + b^ 


(4.14b) 


a/3 = 


-bS 


(4.14c) 


ba = 


-d0 


(4.14d) 


ca = 


—d'j 


(4.14e) 


— 


~cS 


(4.14f) 



Here is a roadmap to the solution: multiply equation (|4.14cp by (|4.14fj) and divide 
by (|4.14d[) and (|4.14e[) to find an expression that can be integrated for 

d = Xivcb/ (va) 
7 = [iiva5/(v/3) 

Substituting these values for d and 7 into the ratio of equations (|4.14d|) and 
(|4.14fj) shows that v — Ai/.iiv, which indicates that v must be separable into a 
product such as v — A2(Z)M2(to), where Ai = A2/A2 and jii = \i<2,jp>i- With the 
above values included and no further integration required, equation (|4.14cp is used 
to find 6 while (|4.14dj) offers c. In summary: 

d = -ba/(3 (4.15a) 

7 = -da/b (4.15b) 

6 = -a/3/6 (4.15c) 

c = -aa/{3 (4.15d) 

The key feature is that, with (|4.15[) . equations (|4.14a[) and (|4.14b[) are automati- 
cally satisfied, furnishing us with no more constraints on the remaining lattice terms 
a, 6, a and (3. Therefore, equations (|4.15|) are the evolution equations, but these 
are not uniquely determined and there is freedom enough to write any equation 
at all into this set, including equations that are known to be not integrable. The 
undeniable conclusion is that any Lax pair of this type is false because one would 
expect that the information to be gleaned about the solution to any evolution equa- 
tion associated with this Lax pair must be underdetermined like equation (|4.15|) . 
Note further that the level of freedom left in the system is exactly that which can 
be removed by gauge transformations, so, in essence, this systems compatibility 
is simply determined by the values of the parameter functions and no evolution 
equation exists. 

4.4. A special case. Item 15 from table [7] is a special case that sees all lattice 
terms linked in all four entries of the compatibility condition. Thus, the compati- 
bility condition supplies only nonlinear equations, none of which can be explicitly 
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integrated, and so no parameter functions present themselves as constants of inte- 
gration. This is an unusual situation but it is not necessary to solve the system 
because of the following arguments regarding the spectral dependence. 

To form the links that define this Lax pair we require the following conditions, 
or an equivalent set, on the spectral terms. 

A = BY/ A E = A 2 /r 

C = BT 2 /A 2 A = A (4.16) 

D = BT/A 

Equation (|4.17p shows the general form of the Lax pair that possesses the links in 
question is 

L - W ( £ »/f> ) (4.17a) 

where Fi(n) = T/A. The prefactors in equation (|4. 17|) are obsolete because they 
cancel in the compatibility condition. This leaves the Lax pair with the same 
dependence on just one spectral term in both the L and M matrices. As such, 
a gauge transformation can be used to completely remove the dependence on the 
spectral variable from the linear problem, implying that equation (|4.17p is actually 
not a Lax pair at all. 

5. CONCLUSION 

All 2x2 Lax pairs, with one separable term in the four entries of each matrix, 
have been considered through the various combinations of terms possible in their 
compatibility conditions. It has been shown that the only non-trivial evolution 
equations that can be supported are the higher order generalizations of the LMKdV 
and LSG equations, with the possible exception of over-determined systems that 
may yet be consistent, see section |4~2H The results of the present work support the 
postulated connection between Lax pairs and singularity confinement [5]. 

There is an important question in where the present work sits in relation to 
studies in multidimensional consistency, or consistency around a cube (CAC), that 
provide a method of searching for integrable partial difference equations where a 
Lax pair can be derived as a bi-product of the procedure [^0, 25 . The present work 
adds to multidimensional consistency studies by removing the restriction that the 
Lax matrices must be symmetric, a consequence of using the same Q equation in all 
three directions of the cube [14] . and also leads to non-autonomous equations, where 
multidimensional consistency studies have only considered autonomous systems. 
Non-autonomy is vital to reductions of the type used in [TB] and [TH] , that lead to 
nonlinear ordinary difference equations and Lax pairs for them. In addition, while 
it has been shown that (CAC) ensures the existence of a Lax pair, the converse is 
not necessarily true, which in itself indicates a need for the present completeness 
study. 

Future studies using the techniques explained here should investigate Lax pairs 
with non-separable terms or with more terms in each matrix entry, and possibly 
aim for a general algorithm for dealing with an arbitrary number of terms in each 
entry of the Lax matrices. Other types, such as differential difference Lax pairs, 
should be considered, as should purely continuous Lax pairs where there is already 
a vast body of knowledge with which to compare results. Further, this technique is 
easily adaptable to isomonodromy Lax pairs. 
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There is some conflict about whether the parameters of the Q4 equation, in the 
ABS scheme, must lie on elliptic curves or not [T? l ITS ! fT§], The Lax pair for Q4, 
found via multidimensional consistency, is for a version of the equation where the 
autonomous parameters are restricted to elliptic curves [20j . An exploration into 
the links present in that Lax pair might resolve the issue by providing a Lax pair 
for Q4 with non-autonomous, and possibly free, parameters. 
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